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Abstract 

We consider the currents formed by a heavy and a light quark within Quantum 
Chromodynamics and compute the matching to Heavy Quark Effective Theory to 
three-loop accuracy. As an application we obtain the third-order perturbative cor- 
rections to ratios of 5-meson decay constants. 



1 Introduction 



Quite often there are problems within Quantum Chromodynamics (QCD) 
involving a single heavy quark with momentum 

p = mv + k , (1.1) 

where m is the on-shell heavy-quark mass and v 2 = 1. In situations when 
the characteristic residual momentum is small (|A; M | <C m) and light quarks 
and gluons have small momenta (|/cf| <C m) it is possible to use a simpler 
effective field theory — Heavy Quark Effective Theory (HQET) (see, e.g., the 
review |T] and the textbooks [2|3] ). Its Lagrangian is a series in 1/m. QCD 
operators are also given by series in 1/m in terms of HQET operators. Here 
we shall consider MS renormalized heavy-light QCD quark currents 

= ZJ 1 (h)jq , jo = qo^Qo , (1.2) 
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where T is a Dirac matrix, and the index means bare quantities. They can 
be expressed via operators in HQET 

j(fi) = Cr(A0KA0 + 2^E^(aOO<(aO + O (^) , (1.3) 

where 

j(fi) = Zj x (ji)% , Jo = qo^Qvo , (1.4) 

<5t,o is the bare HQET heavy-quark field satisfying $Q v o = Q v o, and Oj are 
dimension-4 HQET operators with appropriate quantum numbers. We shall 
not discuss 1/m corrections in this paper; our main subject is the matching 
coefficients Cr (/•«)■ 

The coefficients Cr(/^) have been calculated at one-loop order in the pioneer- 
ing paper 14] . At two loops, they were calculated in Ref. [5], and corrected 
in Ref. [6j|_j All-order results in the large- (3q limit were obtained in [S] (see 
also [7j). Asymptotics of perturbative series were investigated in a model- 
independent way in Ref. [8]. 

In the present paper we calculate the matching coefficients Cr up to three 
loops. These coefficients are useful for obtaining matrix elements of QCD 
currents (such as fs) from results of lattice HQET simulations (see, e.g., 
Refs. j9lfT0] ) or HQET sum rules (see, e.g., Refs. [TUT] ). 



2 Matching 

There are eight Dirac structures giving non-vanishing quark currents in four 
dimensions: 

r = i, 7 «, rj, 

7?7?, JMt, 7?7M, 7W71V, (2.1) 

where 7™ = 7° — ■f)v a 1 and square brackets mean antisymmetrization. The 
last four of them can be obtained from the first four by multiplying by the 
't Hooft-Veltman 7^. We are concerned with flavour non-singlet currents 
only, therefore, we may also use the anticommuting 7^° (there is no anomaly). 
The currents renormalized at a scale \x with different prescriptions for 75 are 
related by [12] 

1 The results incorporating this correction are given by formulae (5.65-5.68) and 
Table 5.1 in |3j. Note a misprint in this table: in the row for T = 7 1 , the term with 
C F in the square bracket, -1453/48, should be +1453/48. 



2 



(>75 AC Q) M = Z P (ji) (qi^O)^ , 
(g 75 AC 7 a Q) M = Z A (jt) (g 75 HV 7 Q Q)^ , 
(^fW 1 ^ = 3t(aO (V^W 1 ^ , (2.2) 

where the finite renormalization constants Zp^a,t can be reconstructed from 
the differences of the anomalous dimensions of the currents. Multiplying T 
by 7^ c does not change the anomalous dimension. In the case of T = ^ a 'y l3 \ 
multiplying it by jf y just permutes its components, and also does not change 
the anomalous dimension, therefore, 

Z T (jt) = 1 ; (2.3) 

Zp,a{p) are known up to three loops [12J. 

The anomalous dimension of the HQET current (11. 4p does not depend on the 
Dirac structure T. Therefore, there are no factors similar to Zp jA in HQET. 
Multiplying T by 7^° does not change the matching coefficient. Therefore, 
the matching coefficients for the currents in the second row of (12.11) are not 
independent: they can be obtained from those for the first row. In the v rest 
frame, where f = 7 , we have 



Z P (fi) 
Za{h) 



Ca C (//) CM 



C HV (H) C 7 7 l 7 2 7 3 
'5 

C 7 AC 7 o(/i) C 7 o(/i) 



Z T (fi) 



C 7 HV 7 o(/i) C 7 l 7 2 7 3(/i) 

C-yAcy^) C^(fi) 

C 7 HV 7 3(/i) C 7 7 l 7 2(/i) ' 

C 7 AC 7 o 7 i(/^) _ C 7 o 7 i(/^) 

C^HV^l {fJ>) C 7 2 7 3(/i) 

= 1 . (2.4) 



- 7^ • 7^ \r~/ -^Yl 

C 7 AC 7 2 7 3(^) _ C 7 2 7 3(/i) 



C 7 HV 7 2 7 3 (/^) C 7 7 l(/i) 

In particular, two matching coefficients are equal: 

C 7x?j (/V) = C^a^ifi). (2.5) 

In the following we shall consider only the matching coefficients for the first 
four Dirac structures in (12.11) . 

In order to find the matching coefficients CV(/i), we equate on-shell matrix 
elements of the left- and right-hand side of Eq. (11.31) . They are obtained by 
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considering transitions of the heavy quark with momentum p = mv + k (11. ip 
to the light quark with momentum k q : 

(q(k q )\j^)\Q(mv + k)) = CM(q(k q )\j(ri\Qv(k)} + O (-, M . (2.6) 

V m ml 



The on-shell matrix elements are! 2 . 



(q(K)\j(v)\Q(p)) = u q (k q )T( P ,k q )u(p) z^zf z]i\ 

(q(k q )\m\Qv(k)) = u q (k q )f(k,k q )u v (k)Z~\ f x)Z 1 J 2 Z 1 q / 2 , (2.7) 

where T{p, k q ) and T(k, k q ) are the bare vertex functions, Zq and Z q are the 
on-shell wave-function renormalization constants of the heavy and the light 
quark in QCD, Zq is the on-shell wave-function renormalization constant of 
the HQET quark field Q v , and Z q differs from Z q because there are no Q loops 
in HQET. The difference between u(mv + k) and u v (k) is of order k/m, and 
can be neglected. It is most convenient to use k = k q = 0, then the 0{\jm) 
term is absent. The QCD vertex has two Dirac structures: 

r(mv,0) —T-(A + Bf) . 

This leads to 

u(p)T(mv, 0)u(mv ) = T(mv, 0) u(0)Tu(mv) with T{mv, 0) = A + B . 

The scalar vertex function T(mv, 0) can be obtained by multiplying the dia- 
grams by a projector and taking the trace. For the first two Dirac structures 
in ( 12. ip . the projector (^+1) can be used; for the next two structures, (^+l)j a . 
The HQET vertex has just one Dirac structure. Therefore, 



T(mv,0)Zr\^Z Q /2 Z^ 
YMZf^Z^Zl 1 



C r(v) - ~,„ '- ^1/2^1% ■ ( 2 - 



Here m is the on-shell mass of the heavy quark (because the external heavy 
quark with p 2 = m 2 should be on its mass shell). Therefore, mass-counterterm 
vertices have to be taken into account on all Q lines. 

If all flavours except Q are massless, all loop corrections to T(0,0), Zq, and 
Z q contain no scale and hence vanish: 17(0,0) = 1, Zq = 1, Z q = 1. If there 
is another massive flavour (charm in 6-quark HQET), this is no longer so. 
The vertex T(0, 0) and Zq have been calculated up to three loops in Ref. [13]. 



2 Both on-shell matrix elements of the renormalized currents are UV-finite. Both 
contain IR divergences, which are the same on the left- and right-hand sides of I 
yielding a finite Cr(/i). 
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The massless-quark on-shell wave-function renormalization constant in the 
presence of a massive flavour has been calculated up to three loops in Ref. [H] 
(an explicit expression can be found in Appendix A of Ref. [13] )• The three- 
loop renormalization of the HQET current Zj has been calculated in Ref. [12] . 

If all flavours except Q are massless, T(mv,0), Zq, and Z q contain a single 
scale m. The on-shell heavy-quark wave-function renormalization constant Zq 
has been calculated at three loops in [IB] (and confirmed in Ref. [T7\); Z q has 
been found in Ref. [T4] . The vertex Y{mv, 0) is the subject of the present pa- 
per. If there is another flavour with a non-zero mass m c < m, there are two 
scales, and calculations become more difficult. The renormalization constant 
Zq has been calculated in this case, up to three loops, in Ref. [18] (the master 
integrals appearing in this case are discussed in Ref. [IS]). The vertex Y{mv, 0) 
and Z q with two masses are considered in Sect. E] and 0J The three- loop renor- 
malization of the QCD currents with all Dirac structures F has been obtained 
in Ref. [20]. 

The QCD quantities T(mv,0), Zq, and Z q in the numerator of ( 12. 8p are cal- 
culated in terms of a° s = g\j (47r) 1_e , the bare coupling of the n/-flavour QCD. 
If there is another massive flavour, say charm, they also contain its bare mass 
m c0 ; we re-express it via the on-shell mass m c . These quantities do not involve 
/i. The MS renormalization constant Zj 1 ^) is expressed in terms of ai nf \fx). 
The HQET quantities T(0, 0), Zq, and Z q in the denominator of (12. 8B are cal- 
culated via a°J and m' c0 , the bare coupling and the bare c-quark mass in the 
^/-flavour QCE0 (we re-express m' c0 via the on-shell mass m c , which is the 
same in both theories). The MS renormalization constant Z~ l (ii) is expressed 

in terms of a s 1 (/i). To combine all these quantities in d 2 . 8 j) . we re-express 

them via the coupling a s f (/x) using the decoupling relation [14] (an explicit 

expression for ai n/ ^(/i) via ai f \fi), including the necessary terms with posi- 
tive powers of e, is given in Eq. (12) of Ref. [21]). 

There exists an exact relation [5] between the matching coefficients C\ (/i) and 
Cj(ji). Namely, the renormalized vector and scalar currents are related by 

id a j a = m(fi)j(ii), (2.9) 

where m(/i) is the MS mass of the heavy quark Q. Taking the on-shell matrix 
element of this equality between the heavy quark with p = mv and the light 
quark with k q = and re-expressing both QCD matrix elements via the matrix 

3 To keep track of flavours in quark loops, we introduce the number of heavy 
flavours with mass m and the number n m of massive flavours with mass m c , so that 
rif = ni + n m + rih and n'j = nj + n m . In reality, rih = 1 and n m = 1. If the c quark 
is considered massless, we can also include it in the number of light flavours n\, and 
set n m = 0. 
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element of the HQET current with T = 1, we obtain 

mCj(ii) = m(/i)Ci(/i) . (2.10) 

The ratio m(fi)/m has been calculated at three loops in Refs. [22.23J. Com- 
paring C^(/z)/Ci(/i) with the analytical result [23] for m(/i)/m provides a 
strong check of our calculations. For m c ^ 0, this ratio has been calculated in 
Ref. [IE]. 



3 Bare vertex functions 



The calculation of T(mv, 0) at m c = is a single-scale problem. The calcula- 
tion is almost completely automated and similar to Refs. [TTf2T] . The Feynman 
diagrams were generated with QGRAF [21] and classified into various topolo- 
gies with the help of q2e and exp [25]. Some sample diagrams are shown in 
Fig. [1] (the first four diagrams). Scalar Feynman integrals were reduced to 
master integrals using integration- by-parts identities [25]. This was done in 
two independent ways: using the Form [27] package SHELL3 [T5] and the C++ 
program Crusher [28] implementing the Laporta algorithm [29]. The master 
integrals for the case of a massless charm quark are known from Ref. [16J (see 
also comments in Ref. [UJ). As an independent check the automatic setup 
described in Ref. [30] has been used and the bare vertex functions have been 
checked numerically with the help of FIESTA [31] . The results for the bare 
vertex functions with all four Dirac structures T can be found on the web 
page [32]. 







Fig. 1. Sample Feynman diagrams for the matching coefficient. Thick, thin and 
dashed lines denote heavy, light and massless quarks, respectively; curly lines are 
gluons. The square denotes the vertex of the heavy-light current. 

Next we study the influence of a non-zero c-quark mass, m c < m, on the 
6-quark matching coefficients Cr- Then c-quark loops exist both in the full 
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QCD and in the 6-quark HQET. The full-theory quantities in the numerator 
of (12. 8p depend on the ratio of the on-shell quark masses, 

x = — . (3.1) 
m 

Some sample diagrams contributing to T(mv, 0) are shown in Fig. [1] (all dia- 
grams except the first four ones depend on x). The HQET quantities in the 
denominator of (I2.8p contain a single scale m c . 

From the technical point of view the calculation is similar to Refs. [T8"f3"3"] . 
We have used Crusher [28J for the reduction. The master integrals are known 
from Refs. [HUH], see also Ref. [31]. Most of the needed terms of their e ex- 
pansions are known analytically, in terms of Harmonic Polylogarithms [351136] 
(HPLs) of x (the status of these expansions is summarized in Tables 1-4 of 
Ref. |19j). From the requirement of cancellation of l/e poles in the match- 
ing coefficients (Sect. E]) we were able to find exact analytical expressions (in 
terms of HPLs) for the O(e ) terms of the master integrals 5.2 and 5.2a (Fig. 8 
in [19]); formerly, they were known only as truncated series in x. The corre- 
sponding entry in Table 3 of Ref. [19] needs updating. We do not present these 
long expressions here, they can be found at 



We have checked the m c dependent results by taking the limit i 4 and 
reproducing the rti part of the results given in Sect. Another check is taking 
the limit x — > 1. If we set = and re-express the renormalized matching 
coefficients Cr via a^ ni \ the results of Sect. [5] with the substitution n m — > 
are reproduced. 



4 Wave-function renormalization of massless quarks 



The last ingredient of Eq. (12. 8p which we had to calculate is the on-shell 
wave-function renormalization constant Z q of a massless quark. The result for 
m c = can be extracted from Ref. [14] . however, we have performed an inde- 
pendent calculation. The m c -dependent part is a new result. The calculation 
is similar to that of Ref. [18], where Zq has been calculated. The integrals 
contributing to Z q reduce to tadpoles when the mass of the incoming particle 
is set to zero. There is only one new type of diagram which is shown in Fig. [21 
all others can be reduced to known results. 

For completeness we give here the result for an arbitrary gauge group. It reads 
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Fig. 2. The diagram which gives the coefficient of n m n/j in Z q . For the meaning of 
the lines see Fig. [TJ 



1 + C F T F (^r(e)^ (n h m~ ie + n m m^) 
16 V 6 36 v ' 



7T 



+ C F Tp — r(e) n ft Z 1 (n ft )m £ + n m Zi(n m )m, 



-6e 



+ T F n h n m (mm c ) 3e Z 2 (m c /m) ) + 0(a£) 



(4.1) 



where a° has the dimensionality m 2e , and the single-scale contributions are [T4TT3] 



Zi(n) 



Z L 
Zh 



CVZp + C^Z^ + T F niZ L + T F nZ H 



e 

96 
1 



443 \ 



12 



e 



1 10 227 2 . , . 

s 1 + — e + e 2 - 16C 3 

192 1 3 9 1 



1879 



e 

72 

e 

36 



l_3 £ + ^ £ 2. 
3 

5 337 
1 e H £ 

6 36 

, 5 151 

1 e H f 

6 36 



+ C»(£ 4 ) 



(4.2) 



where £ = 1 — a , ao is the bare gauge- fixing parameter The two-scale 
contribution (Fig. [2]) is given by 



Z 2 (x) = Z 2 (x- 1 ) = 2Z H -— In 2 x + 0(e a ) 



(4.3) 



4 The gauge-fixing term in the Lagrangian is — (d^A^) 2 /(2ao), the free gluon 
propagator is -(i/k 2 ) (g^ - ffc^fcj, /k 2 ) . 
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5 Matching coefficients 



In this Section, we present the results for the matching coefficients of the 
different heavy-light currents for the colour group SU(3) (results for a gen- 
eral colour group can be found at |32j). For this purpose we decompose the 
coefficients as follows 



C r (n) = l+ aa (m) 4 1) + r S C?\x) 

IX \ 7T 



C?\x) + O(at) , (5.1) 

C r 2) (x) = Cj? + Cf n h + C£ ni + C™{x)n m , 

^(i) = C? G + C? H n h + C^m + C* H nl + C^ L n hni + C^n] 

+ C^ M {x)n m + C^ M (x)n h n m + C^ M {x)mn m + C^ IM (x)n 




2 

m ) 



where T = 1, ^, 7j_, 7j_^, and = rii + n m is the number of active flavours 
in HQET. Furthermore, we use the abbreviation L = ln(/i 2 /m 2 ) (m is the 
on-shell 6-quark mass). 

We first present the results for m c = 0. The individual contributions read 



1 3 2' 

„ 1843 11 2 1 2l n 17 /527 7 2 \ r 9 r2 

Cf = 1 vr 2 H vr 2 ln2 (3 + vr L L 

1 192 72 18 36 s V144 108 J 16 

1 216 18 36 12 

C[ = -^-L 7T 2 --L + -L 2 , 
1 288 36 72 24 



288 36 72 24 ' 

^ rr 8765231 235073 2 917 2l n 1 2l - 3371 

Ci = 1 7r 2 H 7r 2 In 2 it 2 In 2 2 H 

1 62208 46656 324 81 1296 

4733 „ 50039 , 28975 , 4 , 4 32 

H 7r Cs Cfi In 2 a 4 

3888 116640 2592 s 27 9 

/ 46123 25 2 1 2, o 55 95 A. 



■c 



,3 



+ tt 2 H vr 2 In 2 H (3 TT" 

V 1728 108 36 144 Sd 1944 

/ 605 7 2 \ T 2 !5 r3 
+ 1 vr 2 L H Lr 

V 64 48 7 16 
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3349 30917 , 443 2 , n 1 2l 2rt 27845, 29 2 , 

7T 2 + 7T 2 In 2 + -7T 2 n 2 2 Cs + — 7T Ca 

1944 9720 81 9 5184 S3 96 S3 

19 4 45 1 4 8 
-7T — — (, 5 — - In 2 — -a 4 



2430 32^ 9 3 



432 36 6^ J 144 72 

528353 15553 2 25 2l rt 1 2l 2rt 1591, 

7T 2 vr 2 In 2 + — vr 2 In 2 2 Cs 

46656 17496 324 81 972 s 

3281 4 1 , 4 4 / 7399 11 2 97 , 
+ 116640" + 162 ln 2 + 27° 4 + 1-5184 + 972" ~~ 216 



^809 7 r2 19 r3 



+ l864 ~ 648" ) L " 144 L ' 
4045 1 2 11 35 5 2 1 3 

11664 135 27 s3 1296 216 108 
353 1 2 2 u 35 „ 5 1 



5832 216 27" 648 108 54 
6457 13 2 7 , 35 r 5 r2 1 r , 



46656 648 108" 2592 432 216 

- 2 --h, 

3 2 ' 

177 5 2 1 2n „ 11. / 79 7 2 \ r 13 r2 

vr 2 vr 2 In 2 Cs + ?r L + — L 

64 72 18 36 V 144 108 / 16 

727 1 2 

7T , 

432 6 

288 36 72 24 ' 

62575 231253 2 517 2l 20 2 , 2 5645 , 

7r 2 7T 2 In 2 + — vr 2 In 2 2 + Cs 

62208 46656 324 81 1296 s 

2089 2. 17347 4 49435, 11, 4rt 44 

+ vr 2 C 3 vr 4 + — In 4 2 + — a 4 

486 58320 2592 s 54 9 




/ 115 121 2 1 21 n 37, 

+ vr 2 H 7T 2 n 2 H (3 

V 54 648 36 48 S3 

Z2257 91 2 \ r2 13 ~ 
+ (l76 + 432* ) £ "¥ £ - 

2051 _ 24583 361 10 _ 45869 53 

96 2430 27 81 5184 s 96 

1 4 85, 10 , , 80 / 727 1 2 \ r 

7T 4 C 5 n 4 2 a 4 + + — vr 2 L , 

20 32 s5 81 27 V 864 12 / ' 

24457 5575 2 19 2l 1 2l 2 3181 379 4 

46656 + 8748" + 324" ^"Sl" ln 2 + 1944 Cs " 116640" 

1,4. 4 / 319 11 2 83 , 

In 4 2 a 4 + + vr 2 + Cs 

162 27 V 5184 972 216 

/ 469 7 2 \ r2 25 , q 
+ -— -— vr 2 L 2 + 



864 648 / 144 
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5857 1 2 11, 
+ 7777^ + T^Cs , 



7776 405 18 
193 29 2 



432 648 

1751 13 , 7 ,i^r , A L 2_±r3 f53 x 
46656 648 108 s 2592 432 216 ' l ' J 

3 2 ' 

14651 125 o 7 „ 7 , / 31 7 



1728 648 

H L , 

16 ' 

133 5 2 

:7T 



445 1 2 5 1 2 
sei + se 7 " + 72 L ~24 L ' 



2 7 2l rt 7 , / 31 7 2 \ 

IT 2 7T 2 In 2 C3 + 7T 

54 36 s V 144 108 / 



5046967 361033 2 1745 2l n 124 2n 2 ^ 3929, 
-7T 2 - — — vr 2 In 2 + — — vr 2 In 2 2 + — — — C 3 



62208 46656 324 243 1296 

3463 2. 461 4 43835, 163, 4rt 652 

+ tt 2 C 3 7T 4 C5 + In 4 2 + a 4 

972 s 3888 2592 s 486 81 

/301 53 2 T 2 , n 103, 95 4 

+ 7T 2 + 7T 2 n 2 + C 3 7T 

V 54 648 108 144 s3 1944 

fl945 , 91 _ 2 A r2 13 r3 



HW + 432^ 



4133 3385 2 2069 2n 26 2l 2 10445, 

w " «e " + W ln 2 + 2i3" ln 2 " tim 6 

288 14580 96 243 81 

/ 133 5 2 \ r 

+ 1 71 L , 

V 288 108 J ' 

455461 4937 2 169 2l 7 2l 2 5173, 

46666 + 574" + 972" to 2 " 243* ^ 2 + 1944 6 

_ J!»L^ _ X ln < 2 _ 2? a4 + ( _H1 + + ii Cs 

116640 486 81 V 5184 972 216 

+ r864 ~ 648 7r J L + 144 L ' 
3641 11 2 I 7 



+ 77777^ + 77C3 , 



7776 1215 54 
2545 127 2 

-7T 2 



3888 1944 

7993 7 2 7 , 35 r 5 r2 1 r3 

vr 2 C3 H i H L 2 L , (5.4) 

46656 216 108 s 2592 432 216 ' K J 

*-- 5 -L, 
3 6 ' 

20749 7 2 1 2l 5 , / 329 7 2 
— vr 2 - -vr 2 ln 2 - — C 3 + -— - — 7T 2 



1728 24 6 36' V 144 108 



215 r2 
+ 144 L ' 
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^ 648 54 108 36 ' 
7± ^ 2592 108 216 72 

GG 21556403 488167 2 757 2 142 2 8357 

" " 186624 "46656"* ~ 108* 1112 + 243^ ln 2 + 3888 Cs 

16153 2> 2447 4 15925 , 113, 4 904 

H tt 2 C 3 vr 4 (5 H In 2 H a 4 

3888 s 23328 864 s 243 81 

/7871 7 2 5 2l 61. 95 4 

+ 1 7T 2 H vr 2 ln 2 H (3 ^ 

V1728 108 36 48 1944 

/22177 301 2 \ r2 4085 r3 



V 1728 1296 / 1296 



, rH 125005 268333 2 301 2l „ 1 2l - 22469 

T% H . = vr 2 H vr 2 ln 2 + -vr 2 In 2 2 < 

7± ^ 5832 29160 27 9 5184 



59 2 73 4 25. 1, 4 8 

H 7T< 3 7r 4 C 5 In 4 2 a 4 

288 2430 32 s 9 3 



3888 324 18" V 1296 



^ rL 211705 28133 2 71 2l n 1 2l 2n 3347^ 

Cf? j. = 1 vr 2 H vr 2 In 2 vr 2 n 2 2 H (3 

^ 15552 17496 324 27 972 S3 

4183 4 1 4 „ 4 / 4091 13 Q 143 



vr 4 -— ln 4 2--a 4 + — — - — n 2 + — ( 3 )L 



116640 54 9 V 15552 972 216 

/ 3845 7 2 \ r2 5 r3 

+ vr 2 L H Lr , 

V 2592 648 J 16 

H ^ 21281 1 2 31 1 13 2 1 r3 

C ^ " ~34992 + 81 7r + 81 C3 + 144 L + 648 L ~ 324 ' 

a- = -l^ + iI. 2 + A C3 + l L + iiL L2 _^ L 3 ; 

7±? * 17496 216 81 s 72 324 162 ' 

^ LL 29309 89 2 35 /■ 53 r 41 t2 5 r3 ,r rN 

CX L , = vr 2 Cs + ^ + L 2 L 3 , (5.5) 

^ 139968 1944 324 s 2592 1296 648 v ; 

where a 4 = Li 4 (~\. The two-loop results, as well as the coefficients Cf L , are 
known from Ref. all remaining three-loop results are new. 

It is instructive to re-write the matching coefficients for m c = and fi = m 
via the leading /3-function coefficient in HQET, /5q = 11 — |n^. We obtain 



Cj 2j = 0.91/3^ + 1.09 = 7.55 + 1.09 = 8.64 , 
Cf ] = -0.66/3^ + 3.06 = -5.47 + 3.06 = -2.41 , 
Cfl = -1.18/3; + 1.53 = -9.87 + 1.53 = -8.34 , 

CjVf = -1.70/3; + 2.42 = -14.13 + 2.42 = -11.70 (5.6) 
at two loops, and the following results at three loops: 
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C\ 6) = 0.93#, 2 + 9.04/% - 38.16 = 

Cf = -0.54/3^-1.29^ + 29.74 

= -l.28/3' 2 - 5.56/3q + 45.34 

= -l.78{3' 2 - 7.63/% + 63.22 

= -123.96. 



1.74 + 75.34- 38.16 = 101.92, 
-37.25 - 10.72 + 29.74 = -18.23 , 
-88.92 - 46.34 + 45.34 = -89.92 , 
-123.61 - 63.57 + 63.22 

(5.7) 



A method to estimate higher loop contributions called naive nonabelianization 
has been formulated in Ref. [5j. It is based on the fact that each polynomial 
in rif can be re- written as a polynomial in fl . Usually it is relatively simple to 
calculate the term with the highest power of rif. This means that we know the 
coefficient of the leading power of (3 . Neglecting sub leading powers of (3 we 
obtain an estimate of the full result. The two-loop corrections to the matching 
coefficients (15.61) were among the examples confirming naive nonabelianiza- 
tion [5]. At three loops we see that this prescription reproduces the correct 
signs and roughly the correct magnitude of the full results. In all cases we ob- 
serve a compensation between the O(/3' ) and 0(1) terms. In the case of 

(3) 

and this compensation is almost complete, and naive nonabelianization 
works surprisingly well. 



Since the expressions for the charm-mass dependence are quite involved and 
not completely expressed in terms of HPLs of x (the 0(e) terms of the master 
integrals 5.2, 5.2 known analytically only as truncated series 

in x, see Table 3 in Ref. [19]), we refrain from listing them here. Instead, we 
present an expansion of our results to the second order in x. The results in 
terms of the master integrals and expansions to higher orders can be obtained 
from the web page [32J. 

Our results for the m c -dependent coefficients read 
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6 Meson matrix elements 



We are now in the position to apply our results to the matrix elements between 
a B or B* meson with momentum p and the vacuum. They are defined through 



(0| (qis°Q) J£> = -im B fE(»), 

(0\qn£ C Q\B) = if B p a , 
(Olq^QlB*) = im B *f B *e a , 

(0| (q^Q)^ \B*) = fl(fi)(p a e^-p^e a ) , (6.1) 

where e a is the B* polarization vector. The corresponding HQET matrix ele- 
ments (at m oo) in the v rest frame are 



q>y£ c Q v ) \B{h)) m = -iF(fi), 
(0|(^)jS*(^)) nr = zF(/i)e, (6.2) 

where the single-meson states are normalized by the non-relativistic condition 

jB(k')\B(k)} nr = (27r) 3 5(k'-k). 

We also remind the reader that qT$Q v = qTQ v , so that there are only two 
currents. These two matrix elements are characterized by a single hadronic 
parameter F{p) due to the heavy-quark spin symmetry. From Eq. (12.91) we 
have [5] 

/IW _ m B 



Ib rn(fj,) 



(6.3) 



Here m B = m + A + (9(AQ CD /m) where A is the residual energy of the ground- 
state B meson in the limit m — > oo. Neglecting 1/m corrections, we see that 
this equation coincides with (I2.10p . We have checked that our results agree 
with the known formulas for m(/i)/m at m c = [23] and with m c corrections 
taken into account [Tg] . 

Using the expressions of Section we find the following results for the ratios 
of decay constants. Again, we present an expansion to the second order in x 
for the charm mass dependent terms. For the numerical evaluations we use 
an expansion to the eighth order in x and the values a^^rtib) = 0.2163 and 
x = 0.3. 
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(6.4) 



In the second line from the bottom the corrections from tree level, first, second 
and third order in a s are given separately. Also the contributions stemming 
from the finite charm mass are separated (the second number in the paren- 
theses). In the first part of the last line, the m c correction is also separated. 
Power corrections O (Aqcd/ nib) are discussed in Refs. [3~Tf8] and amount to 
several per cent. 

For the second ratio we obtain for /i = m b 
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(6.5) 



The coefficients of these perturbative series (at m c = and /i = m b ) can be 
re- written via /3' : 



( /b* 
' fl*{m b ) \ 

, fB* ) 
\fB. 

7b* ( m b) 



(2) 



13* 



(2) 



(3) 



(3) 



-0.53/3q - 1.97 = -4.40 - 1.97 = -6.37. 

-0.51$ + 0.89 = -4.26 + 0.89 = -3.37 , 

-0.74/?; 2 - 5.06/?; + 16.33 

-51.67 - 42.21 + 16.33 = -77.55 , 

-0.50/3q 2 - 2.75$ + 19.07 

-34.69 - 22.91 + 19.07 = -38.53 . 



(6.6) 



Again, naive nonabelianization [5] works quite well for these ratios predicting 
the correct sign and order of magnitude. 
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Asymptotics of the perturbative coefficients for the matching coefficients at 
a large number of loops I ^> 1 have been investigated in Ref. [8] in a model- 
independent way. The results contain three unknown normalization constants 
^0,1,2 ~ 1- Let us in the following assume that the number of loops I = 3 is 
much larger than one and compare our results with Ref. [Sj. The asymptotics 
of the perturbative coefficients for /b*//b contain N Q and N 2 (see (5.6) in 
Ref. [8]); in the case of m/rh it contains only Nq (see (5.9) in Ref. [8]1 5 1: 



'f_B_ 
Jb ) L 



(n+1) 
5/3 



14 
"27 



O 



50 



-n 



-9/25 r 



o 



N2 
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/ m \(n+l) 



\rhj 



L=-5/3 



(6.7) 



The coefficient of ./V2/./V0 is about 0.08 at n — 2, and it seems reasonable to 
neglect this contribution. Neglecting also 1/n corrections, we obtain [5] 



7b. 

Jb 



(3) 



-5/3 



14 

-— • 56.37 
27 



-29.23. 



(G.i 



Our exact result 



Jb* 
Jb 



(3) 



! = -37.787 (6.9) 

L=-5/3 

agrees with this prediction reasonably well, thus confirming the simple rela- 
tion (5.14) in Ref. [8]. However, 1/n corrections are large and tend to break 
this agreement. It is natural to expect that 1/n 2 (and higher) corrections are 
also substantial at n = 2. 



7 Conclusion 



We have calculated the N 3 LO corrections to the matching coefficients of heavy- 
light currents in HQET. Our result takes into account effects due to the mass 
of the charm quark. Strictly speaking, our results should be used together with 
the N 3 LO /3-function [381139] and anomalous dimensions of both the QCD cur- 
rents and the HQET one. The four-loop anomalous dimensions are known for 
some of the QCD currents (for r = 7j_ the anomalous dimension is exactly 
zero; for T = 1 it is just the anomalous dimension of the MS mass [40IT4T] with 
a minus sign). However, the four- loop anomalous dimension of the HQET cur- 
rent is not known (this anomalous dimension does not depend on the Dirac 

5 Note that for convenience the choice [i = me -5 / 6 has been adopted in Ref. [8]; m 
is the renormalization-group invariant mass (the exact definition used here is given 
in the unnumbered formula after (3.8) in Ref. [8]). 
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structure r)0 The effect of this unknown anomalous dimension cancels in ra- 
tios of S-meson decay constants, fs*/ fs and /b*(aO//b*, discussed in Sect.[6j 

All the previous experience shows that contributions from anomalous dimen- 
sions are numerically much smaller than from matching coefficients. Matching 
coefficients have renormalon singularities at the Borel parameter u — 1/2, this 
is the position closest to the origin out of all possible ones; this means that the 
factorial growth of their perturbative coefficients is fastest among all possible 
variants. On the other hand, it is generally believed that anomalous dimen- 
sions have no renormalon singularities, and their perturbative series have finite 
convergence radii. 



Only fs has been measured experimentally [42] . Our results can be used for 
predicting the B* decay constants. We find that the perturbative series for 
/b* / fn and /J, / converge very slowly at best. The effects due to the charm- 
quark mass are small and of the order of 10 -3 . 

The matching coefficients Cr can be used for extraction of B (and B*) decay 
constants from lattice HQET simulations (see Ref. p3] for recent reviews), or 
from HQET sum rules. 
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